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BEHAVIOR  OF  SEPARATION  CONSTANTS 
FOR  FINITE  GYROMAGNETIC  PLASMAS 


Abstract 

The  type  of  the  functions  appearing  in  the  expression  for 
electromagnetic  fields  in  finite  magnetoplasmas  depends  only 

on  the  properties  of  the  medium  manifested  in  the  separation 

« 

constants.    The    co-j3   plane  of  the  Brillouin  diagram  is 

partitioned  into  a  definite  number  of  zones  characterized  by 
the  type  of  field  solutions.    The  general  quantitative  behavior 
of  the  Separation  constant  for  finite  magnetoplasma  in  a  longi- 
tudinal magnetic  field  on  the  Brillouin  diagram  is  described. 


BEHAVIOR  OF  GYROTROPIC  PLASMA  SEPARATION  CONSTANTS 
AS  FUNCTIONS  IN  THE    o>-|8    PLANE 

Of  the  numerous  papers     concerned  with  the  propagation  of  waves 

through  waveguides  filled  with  gyrotropic  plasma  columns  the  majority 

2   3   4  5 

is  confined  either  to  the  qua  si- static  approximation    '    '    '     or  to  a  dis- 

7 

cussion  of  the  problem  in  general  terms  with  tabulation  of  auxiliary  functions    , 

o 

such  as,  e.g.  ,  the  biradial  functions    and  an  outline  of  the  general 
method  of  approach  for  obtaining  the  solution    '    '    '       without  a 
systematic  and  comprehensive  presentation  and  classification  of  the 
numerous  types  of    a."  -  #   propagation  diagrams  and  field     solutions  such 

as  are  known,  e.g.  ,  for  empty  waveguides  or  waveguides  filled  with  the 

11  12 

dielectrics       and  for  slow  wave  structures.         Only  a  few  representative 

solutions  have  been  published  thus  far      '      '       and  they  did  not  include 

a  complete  systematic  survey  of  all  possible  cases  although  they  were 

certainly  oriented  in  this  direction. 

Despite  a  high  number  of  combinations  that  can  arise  with  the 

several  parameters  that  can  be  varied,  viz.  ,  the  plasma  density  or  the 

corresponding  plasma  frequency    w     ,  the  intensity  of  the  collimating 

P 

magnetic  field  or  the  corresponding  cyclotron  frequency    co       ,  and  the 

cy 

diameter  of  the  waveguide  or  the  corresponding  empty  waveguide  cutoff 
frequency,  it  is  nevertheless  possible  to  establish  a  system  of  classifi- 
cation of  modes  of  propagation. 

It  is  believed  that  one  way  to  make  it  possible  readily  to  obtain  the 
actual  solutions  for  given  conditions  in  waveguides  is  to  prepare 


catalogues  where  quantitites  of  interest  such  as,  e.g.  ,  the  arguments  of 
Bessel  functions,  biradial  functions,  and  the  like  are  tabulated  and 
presented  graphically  on  charts  and  diagrams  as  functions  in  the    w~  $ 
plane  for  various  combinations  of  values  of  the  plasma  frequency, 
cyclotron  frequency,  and  the  cutoff  frequency  of  the  empty  waveguide. 
In  addition  to  these,  actual  solutions  for  combinations  of  the  above 
values  in  reasonable  and  suitably  chosen  steps  would  also  be  catalogued. 

The  objective  of  this  paper  is  to  present  a  study  of  the  separation 
constants,  if  indeed  these  quantities  may  be  called  by  this  name,  which 
in  circular  cylindrical  geometry  appear  as  arguments  of  Bessel  functions 
in  the  expressions  for  the  electromagnetic  fields    '    '    '      ,  viz.  , 


clIn<uir)  +  c2In(VT) 


\  -  2         2  V 

Ul'u2 


'rti'Ar^v'W"!!1 


Hz 2"2  (2) 

U1"U2 


Evidently,  the  separation  constants  determine  the  kind  of  the  Bessel 

functions  in  Eqs.   1  and  2.     If  either    u      or    u      or  both  are  imaginary, 

the  corresponding  Bessel  function  become  modified,     I      instead  of    J   . 

n  n 

If  they  are  complex  conjugates  Beseel  functions  of  complex  arguments 
appear  in  the  above  expressions.    It  is  therefore  believed  that  it  is  more 


useful  to  know  the  behavior  of  the  separation  constants    u      and    u      in 
the    60-/3  plane  than  merely  the  behavior  of  the  separate  components 
of  the  equivalent  dielectric  tensor    c  which  is  sometimes  presented. 

The  constants    u1    and    u?    are  solutions  of  a  biquadratic  equation 
which  is  obtained  solely  from  the  properties  of  the  medium.     In  addition, 
of  course,  they  must  like,  c      and    c„,  also  be  compatible  with  the 

boundary  conditions.    The  medium  or  electronic  equation  was  first 

g 
derived  by  Hahn    and  we  shall  give  it  in  its  original  form: 


[Cf)2^2^2]2=^r^f)2^2-2][Cf)2<-  jW-2)]   « 


The  following  notation  is  being  used  here: 


k    = 

P 

c 

"p  = 

J- 

k      = 

cy 

c 

(4-a) 


(4-b) 


(4-c) 


»cy="o    frH0  (4~d> 


k=    ^  (4-e) 


Symbols  denoting  the  quantities  are  as  follows: 

co      =      frequency  at  which  the  waves  in  waveguide  propagate; 

pQ     =      electronic  volume  charge  density; 

jS       =  propagation  constant  appearing  in  the  expression    exp  j(a>t-jSz); 

c       =  velocity  of  light; 

e/m  =  1.76  x  10       coulomb/Kg,  electron  charge  to  mass  ratio; 

H.     =  axial  static  magnetic  field  strength; 

—  fi 
jlu  =  1.25  x  10       henry/m; 

-12 
E       =      8.86x10         farad/m; 

i 

a       =      radius  of  the  cylindrical  waveguide. 

The  explicit  solution  of  the  biquadratic  equation  is: 


P     cy 


2      2  / 

k     k  /  2 
cyp            /„2  ,2.2    An2    k        „2  ,2,1  ... 
\/(j3   -k  )  +40     — =-  (k  -k  )  \                                       -   (5) 

k"  k  p  J 

cy 
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The  plus  sign  in  front  of  the  square  root  will  be  consistently  used  for  the 

2  2 

calculation  of    u.    and  the  negative  sign  for  the  calculation  of   u„. 

2 
It  is  evident  that  we  may  multiply  the  Equation  (5)  by    a      throughout 

2 

and  then  obtain    u      expressed  in  terms  of   ka,  k  a,  k     a,  fia     which  is 

P         cy       H 

useful  for  representation  on  diagrams  since  these  values  are  then  pure 


numbers . 

2 

The  following  relationship  holds  between  the  function    u     and  the 

2  10 

function    x      used  by  Suhl  and  Walker     : 


2 

2  Ul,2  ,.. 

Xl   2  = 2 2  (6) 

lt2      (ka)Z-(k  a)Z 


It  appears  that  the  reason  for  the  use  of  this  expression  by  Suhl  and 
Walker  was  the  desire  for  unified  treatment  of  both  ferrite  and  plasma 
filled  waveguides.    Apparently  propagation  was  not  expected  for  values 
of  frequency  lower  than  the  value  of  plasma  frequency.    A  numerical 

solution  of  the  equations  of  Suhl  and  Walker  would  have  discovered  the 

2 
plasmaguide  modes  that  were  later  found  by  Trivelpiece  and  Gould   . 

14 
Camus  and  Le  Mezec      use  the  notation: 


T2  =  -^-  (7) 

a 

and 

2 

Z  =  — ^  (8) 

(kpa)2 


When  boundary  conditions  are  introduced  it  is  best  to  normalize  all 
frequencies  to  the  frequency  defined  by: 


WQ  =  c/a  (9) 


This  means  that  quantitites    k  a,  k      ,  ka    can  be  considered  essentially 

P         cy 

as    co   ,  co      /to    normalized  with  respect  to    ton.    Whenever  boundary 
p       cy  U 

conditions  are  introduced  the  radius  of  the  waveguide  appears  in  the 

expression  and  both  plasma  frequency    co     and  cyclotron  frequency    co 

are  specified  by    k  a    and    k     a    respectively,  i.e.,  in  terms  of  the 

P  cy 

normalizing  frequency    con  =  c/a .    The  solution  of  the  boundary  condition 

problem  depends  on  the  absolute  value  of    co      and    co      .    However,  when 

P  cy 

2  2 

we  consider  only  the  quantities    u,    and    u?    as  functions  of  the  frequency 

co   and  wave  number    jS  ,  i.e.  ,  as  functions  in  the    co-j3    plane  it  is 

possible  to  normalize  all  frequencies  to  the  plasma  frequency  (or  alterna- 

tively  to  the  cyclotron  frequency)  and  the  ratio  of  the  cyclotron  frequency 

to  the  plasma  frequency  becomes  the  only  parameter  that  is  left  to  be 

varied.    This  ratio  is  denoted  in  the  following  way: 


co  (k     a) 

M  =  — SL.    =       <2_  (10) 

«p  (kpa) 


Likewise,  it  is  convenient  to  normalize  the  other  quantitites  with  respect 
to  the  plasma  frequency: 


Y_     (ka) 
(kpa) 

CO 

(k  a) 
.  P 

CO 

p 

(11) 


(12) 


2  2 

In  this  normalization  the  Hahn's  equation  for    u.    and    u.    reads: 


2 
u,    „  ,  ,,2 


u?  .--^ 


i r{-[2Cf2-1-MV  -^]0<2-Y2)-2(Y2-l) 


1,2     (k  a)2      2(Y2-1-M*V    -  Y 


Y  M 


(13) 


2  2 

Note  that    U1    and    U9    are  multivalued  at  the  origin    (X=0,Y=0).     Its 

value  depends  on  the  path  in  the    Y-X    plane  along  which  we  approach  the 

2 
origin'.    It  is  seen  without  difficulty  that  the  values  of    U  are  real 

1  |  6 

for  values  of  X   and   Y   exterior  to  the  region  in  which  the  expression 
under  the  square  root  becomes  negative ,  viz . , 


0  <  Y   <    1  (14) 


^[/l+M2-Y2      -   /l-Y2    ]<X<-^[Vl+M2-Y2  +    n/i-Y2    ]  (15) 


2  2 

Within  this  region  the  values  of    U.    and    U„    are  complex  conjugates. 

The  boundary  of  the  complex  region  is  a  closed  curve  touching  the  line 


1+M                       / 1  +  M 
Y  =  1    at   X  =  1    and  the  line   X  =   / —-    at   Y  =  /- ~-      .As    M  is 

M  V  2+M 

increased  the  complex  region  becomes  narrower.    This  boundary  is  also 

shown  on  Figures    1  to  4  and  6  to  15.    On  the  above  boundary  curve 
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Fig»   1    Contour  curves  for  constant  values  of    Re(U  )    in  the  complex 
region  of  the    ts  -  jS  plane  at  the  value  of    M  =  0o8o 


Fig,  2    Contour  curves  for  constant  values  of    Im(U  )    in  the  complex 
region  of  the    w  =  0   plane  at  the  value  of    M  =  0„8„ 


10 


Fig.  3    Contour  curves  for  constant  values  of    Re(U  )    in  the  complex 
region  of  the    w    -    (5   plane  at  the  value  of    M  =  1.6. 


11 


2 
Fig.  4    Contour  curves  for  constant  values  of    Im(U  )    in  the   complex 

region  of  the    05  -  0  plane  at  the  value  of    M  =  1.6. 


12 


2 

(Eqs.   14  and  15)  the  expression  under  the  square  root  vanishes  and    U,  = 

2  2         2 

U„.    The  values  of    U.  =  U9    on  the  boundary  curve  are: 


(16) 


on 

X 

Y2 
M2 

=  3L 

M 

■'[/l  +  M2-Y2 

i-r  m2 

L2Y2(1+M2-Yl 

-] 

/> 

-2        ] 

U2=2 

/(M2- 

fl-Y2)(l-Y2) 

(17) 


and  on 


x-£  {J^F^      +  J^] 


(18) 


U^IJY2-!^ — -   -1,1     y(M2+l-Y2)(l-Y2)    }  (19) 

MZ    ^  L  2Y   (1+M  -Y  )        J  '  J 


If  the  origin    00=0,  Y=0)  is  approached  along  the  curve  defined  in  Eq.   16 

2 

the  value  of    U     at  zero  is  given  by 

U2  =  -  — - (20) 


/ 


1+M2 


if  on  the  other  hand  the  origin  is  approached  along  the  curve  defined  in 

2 
Eq.   18  the  value  of    U      at  zero  is  given  by 


U2=+  — - (21) 


/• 


1+M 


13 


The  equations  of  contour  lines    Re(u  )  =  const  and    Im(IT )  =  const  within 

2 

the  region  defined  by  Eq.  14  and  15  are  readily  found.    For  Re(U  )  =  constant 


X=     /y2+2      (1-yV(1+M2-yV(U2)  (22) 

[    M.  -2(m2+1-Y2)  ] 


Y 


2 
for    Im(U  )  =  constant 


X=  X-  y(M2+l-Y2)+(l-Y2)±2    y(l-Y2)(l+M2-Y2)-(M2+l-Y2)2(Im(U2)  )2 


The  curves  given  by  Eq.  22  are  single  valued  functions  of   Y   while  the 
curves  given  by  Eq.  23  are  double  valued;  they  are  closed  curves 
similar  to  the  boundary  of  the  region.    Figures  1  to  4  show  representative 
contour  diagrams  inside  the  complex  region.     It  is  further  of  interest  to 

determine  whether  and  in  which  points  on  the  boundary  curve  the  value 

2  2 

of    U     vanishes.    Evidently    U  =0   at  Y=l    and  X=l ,  other  points  are 

found  by  setting  the  expressions  of   Eq.   17  or  19  equal  to  zero,   and 

2 

arrange  the  terms  in  powers  of   Y  .    Thus  we  obtain: 


Y6-(2+M2)Y4+(l+M2)Y2  -  -~  =  0  (24) 


The  roots  of  this  equation  are  real  as  would  become  apparent  if  Eq.  24  were 
solved;  it  is  possible  to  obtain  them  by  solving  them  numerically.    However, 


14 


a  more  elegant  way  to  study  its  roots  is  available.    It  will  be  recalled  that 

2  2 

either    U      or    U9    can  vanish  on  the  following  two  curves: 


X=       hf -T7  (25) 


Y2- 

l 

Y 

'*- 

1 

Y 

X=    VY    -— M  (26) 

1+T 

These  two  curves  lie   outside  of  the  complex  region.    Equation  25  re- 
presents a  curve  with  a  discontinuity  at  Y  -  M,  its  upper  branch  inter- 
Mi         /  M^ 
cepts  the   Y  axis    (corresponding  to    co   axis)  at  Y=  +  -r-  +  »/—■ - —    +  1; 

Eq.  26  represents  a  continuous  curve  with     Y    intercept  at 


Y  =  — —  +  J  -j-     +  1.    Values  of    U.    and    U?    on  curves  given  by 
Equations  25  and  26  can  be  readily  evaluated.     On  the  curve  given  by 
Eq.  25  we  have: 


U2  M  r  Y3-MY2-Y+m4  ,       If    Y3-MY2-Y+M/2n2h  (2?) 

1,2    Y^-l-M2)'"         Y-M  ~   V  ^         y-M  ■     ~ 


Care  must  be  exercised  in  taking  the  square  root  in  this  expression. 
We  shall  here  take  the  square  root  as  the  opposite  operation  of  squaring. 
The  expression  under  the  square  root  is  always  positive  but  the  expression 
squared  can  be  either  positive  or  negative.    For  the  sign  of  the  square 
root  we  take  the  sign  of  the  squared  quantity  which  we  must  first  evaluate. 
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Fig.  5    Plot  of  the  parameter    M    as  a  function  of   Y(=o>Va5p)  illustrating 
the  solutions  of  the  cubic  of  Eq„  28  „ 
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The  roots  of  the  cubic  equation 


Y3-MY2-Y+  -j-  =  0  (28) 


may  be  determined  as  follows.    From  Eq.  28  we  have 


__         Y(Y-1)(Y+1) 

M 1 i 

(Y-=4)(Y+^:)  (29) 


The  value  of    M    can  be  readily  plotted  as  a  function  of   Y   as  is  shown 
on     F*ig.  5         .    From    Fig.  5        it  is  seen  that  one  root  of  Eq.  28  de- 
signated as  first  root,  is  always  negative  and  that  the  third  root  is 
always  found  between  the  following  limits: 


M<Y3<-f  +  J(-f)2+l  (30) 


.  If  we  limit  our  interest  to  the  first  quadrant  of  the    to-/S    plane 
(Y  >0,  X  >  0)    and  to  real  values  of  X   these  two  roots  need  not  be 
considered.    The  second  root,  on  the  other  hand,  lies  between  the  values 


°  *  Y2    Vf  <31) 


This  means  that  a  point  exists  on  the  lower  branch  of  the  curve  given  by 

Eq.  25  with  an  ordinate    Y   which  satisfies  the  cubic  of  Eq.  28  and  as  a 

2  2 

consequence  both    U.    and    U?    are  equal  to  zero  in  that  point.    If  both 


17 


2  2 

U.    and    U„    are  equal  to  zero  then  they  are  also  equal  to  each  other  and 

consequently  lie  on  the  curve  given  by  Eq  ,    18,  i.e.  ,  on  the  right 

2  2 

boundary  of  the  region  of  complex  values  of    U,    and    U9 .     Hence  the 

curve  of  Eq.  25  has  one  point  in  common  with  the  right  boundary  of  the 
complex  region  given  by  Eq.   18  and  the  two  curves  are  tangent  to  each 
other  in  that  point.    We  shall  call  this  point  the  lower  tangent  point.     Now 
for  values  of   Y    between  zero  and  the  second  root   Y?    of  Eq.  28  the 
squared  expression  in  Eq.  27  is  negative,  viz.  , 

'  Y3-MY2-Y+  -y-   >     0  (32) 


Y  -  M  <  0  (33) 


if  0  <  Y    <Y2  (34) 

2 
From  where  follows  for    U. 


u?=    M 


2  2 

Y(Y  -1-MZ) 


f  Y  -MY-Y+M^2  ,  r  Y  -MY-Y+M/2  "N  1      n  ,%tA 

i    y-m  C    y^m        ;r°  {6b) 


2 

and  for    U_: 


,3 


u2=    2M(Y  -MY-Y+  M/2) 
2     Y(Y2-1-M2)(Y-M) 


(3  6) 
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When,  on  the  other  hand,  the  value  of   Y    is  greater  than  the  value  of  the 
second  root  the  following  relations  hold: 


3-MY2-Y+-^-  <  0  (37) 


Y 


Y-M  <  0  (38) 


Y2  <    Y    <  M  (39) 


The  squared  quantity  in  Eq.  27  is  then  positive  and  we  have 


U2  =     2M(Y3-MY2-Y+M/2)  (4Q) 

1         Y(Y2-1-M2)(Y-M) 


TT2  M f  Y3-MY~Y+M/2      Y3-MY2-Y+M/2 1      n  .:.. 

u* = vw, "^        s=s     > =0  (41) 

Thus  it  is  seen  that  on  the  curve  given  by  Eq.  35  between  the  origin 
(X=0,  Y=0)  and  the  point  in  which  this  curve  touches  the  boundary  of  the 
complex  region  the  function    XL     is  zero  and  the  function    U7    has  a  value 

given  by  Eq.  40,  on  the  rest  of  the  curve,  beyond  the  lower  tangent  point, 

2 

the  situation  is  reversed  and  the  function    U^    is  zero  while  the  function 

2 

U.     is  given  by  Eq.  40.     Likewise  it  can  be  seen  that  on  the  upper  branch 

2 
of  the  curve  given  by  Eq.  25  the  value  of  the  function    U      is  always  zero 

2 
and  the  value  of  the  function    U      is  given  by  Eq.  40  for  positive  values 

19 


of  Y. 


2  2 

In  the  same  manner  we  obtain  the  values  of    U.    and    UL    on  the 


curve  given  by  Eq.  26. 


TT2  M f     Y3+MY2-Y-M#,     /'Y^j^Y^-Y-M'fc  V\ 

V"  ^^2j   I"  — Y^M ±A YTM J   I  (42) 


To  determine  the  sign  of  the  square  root  it  is  again  necessary  to  study  the 
behavior  of  the  cubic  polynomial  and  the  roots  of  the  equation 


Y3+MY2-Y-  -y-=  0  (43) 


It  is  immediately  obvious  that  this  equation  becomes  identical  to  Eq.  28 
if  the  variable   Y    is  replaced  by    -Y.    The  roots  are  then  as  shown  on 
Fig.      5         only  with  the  positive  and  negative    Y   axis  interchanged. 
From  the  same  considerations  as  before  only  the  positive  root,  designated 
henceforth  as  the  third  root  (of  Eq.  43)  is  of  interest.    The  third  root, 
lies  between  the  values 


P=  <  Y.    <  1  (44) 


This,  again,  means  that  a  point  exists  on  the  curve  given  by  Eq.  26  with 

an  ordiante   Y-    which  satisfies  the  cubic  of  Eq.  43  and  consequently  both 

2  2 

U1    and    U„    are  equal  to  zero  in  that  point  and,  according  to  the  argument 
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_ — ,_-_ p_1 __.    —,,^_,..  TTT" r"  ' 


presented  before,  this  point,  too,  lies  on  the  curve  given  by  Eq.   16,  i.e.  , 

2  2 

on  the  left  boundary  of  the  region  of  complex  values  of    U,    and    U« .    Thus 

the  curve  of  Eq.  26  has  a  point  in  common  with  the  left  boundary  of  the 
complex  region  given  by  Eq.   16,  the  two  curves  are  tangent  to  each 
other  at  that  point.    We  shall  call  this  point  the  upper  tangent  point.    Now 
for  values  of  Y   between  the  ordinate  intercept  of  the  curve  given  by  the 
Eq.  26  and  the  third  root,  viz.  , 


M  +   If  M  ^2 


ir)2  +1     <  Y  <  Y3  (45) 


the  squared  expression  in  Eq.  42  is  negative 


Y3+MY2-Y-  y  <    0  (46) 


2 

From  where  follows  for    U. 


1     Y^-l-M2)    L         Y+M  J 

2 
and  for    U 


TT2  M r      Y3+MY2-Y-M/2     s      Y3+MY2-Y-M/2^  -i 

2=     Y(Y2-!-M2)    t"— ^M C-— YT^ )]-»■  <48> 
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when,  on  the  other  hand,  the  value  of   Y   is  greater  than  the  value  of  the 
third  root  the  following  relations  hold 


•Y3  +  MY2  -  Y  -  -f-   >  0  (49) 


Y>Y3  .  (50) 


2 

The  squared  quantity  in  Eq.  42  is  then  positive  and  we  have  for    U. 


TT2        M                 r      Y3+MY2-Y-M^,  /Y3+MY  -¥-M/2>  1       n  tc,\ 

l=W^Y~™i +C       Y+M  jj=°  (51) 

2 

and  for    U_ 


TT2  _      -2M  f    Y  +MY  -Y-M/2  1  ,„\ 

2 2  2~  Y+M  {     * 

1      Y(Y  -1-Ml     L  Y   M  J 


In  the  same  manner  it  can  be  shown  that  on  the  left  boundary  of  the 

2 

complex  region  given  by  Eq.  16  the  values  of    U     are  negative  for  values 

of   Y   samller  than  the  value  of  the  third  root    (Y<  Y«)    and  positive  for 

values  of   Y   greater  than  the  value  of  the  third  root    (Y_  <  Y  <  1)    and 

2 

that  on  the  right  boundary  given  by  Eq.   18  the  values  of    U    are  positive 

for  values  of   Y   smaller  than  the  value  of  the  second  root  (Y  <  Y?)    and 
positive  the  values  of  Y   greater  than  the  value  of  the  second  root 
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(Y    <  Y   <  1).    Moreover  it  can  be  seen  that  on  the  line    Y=l    the  values 

2  2 

of  functions    U.    and    U„    are  as  follows: 


In  the  interval: 

Y=l, 

0  < 

X 

< 

1 

<- 

0 

<- 

1  - 

x2 

(53) 


(54) 


(55) 

and  in  the  interval: 

Y=1(KX<    •  (56) 


U*    =  1  -X2  (57) 


U22  =  0  (58) 


To  show  that  the  roots  of  Equations  28  and  43  are  also  the  roots  of 

2  2 

Eq.  24  obtained  for  the  condition  that  both    U      and    U9    vanish  we 

multiply  the  cubic  polynomials  of  Equations  28  and  43;  the  resulting 

bicubic  polynomial  is  just  that  of  Eq.  24. 

2  2 

Since  in  the  complex  region  the  functions      IL    and    U-    are 

complex  conjugates  they  can  be  represented  by  their  modulus    p      and 

argument   <p     ,  viz.;      U.  =  p  /  <p;    U„  =  p  /  -cp    .    Note  that  we  have 

9  9 

here    U1    and    U2    instead  of    U^    and    U^.    Contour  curves  of  constant 

23 


modulus  and  argument  of  the  functions    U     and    \J     are  of  greater 

practical  value  than  the  contour  curves  for  constant  real  and  imaginary 

P  ■  1  fi 

part  of    U   i    Tables  of  Bessel  functions  of  complex  arguments       give  the 

argument  of  the  Bessel  functions  in  polar  form. 

The  modulus    p    is  readily  found  from  Eq.   13: 


2 

r  4(Y2-1-M2)2    lL  YZ  J 


+  4  T4    XV     (1.Y2h(x2.y2)2  -  | 


To  find  the  equation  of  a  contour  curve  in  the    Y-X  (<*}  -  /3  )    plane  for 

2 
constant  value  of    n     the  above  equation  is  expanded  in  powers  of   X 

2 

and  a  quadratic  equation  in   X      is  obtained. 


C4X4  +  C2X2  +  CQ  -  pA  =  0  (60) 


The  coefficients    C   ,  C« ,  C.      depend  only  on   Y  and    M   and  are  given 
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by  the  expressions: 


V 


[20f2-l-M2)+  ^V^yV-D^-i-mV  $>4U-y2)2-m4 

4(y  -i-M  ) 


(61) 


C2= 


-2Y2[2(Y2-1-M2)+|^]2+4(Y2-1)[2(Y2-1-M2K  ^^(l-Y2)^ }   (62) 


2  2  2 

4(Y  -1-mV 


^     _  'Y2-l-M2+M2/Y2 

G4  "         2  2 

Y  -1-M 


(63) 


From  Eq.  60  we  obtain  the  explicit  form  of' the  equation  for  the  contour 
curve    of  constant  modulus  -  p   • 


X  = 


2C4    " 


H) 


2  *4  - 

2+      "*T=- 


Cq 


(64) 


Obviously   X    must  be  real  and  positive.     Some  representative  contour 
curves  given  by  the  above  equation  are  shown  on  Figures     6  and    8 . 
The  argument   <p    is  also  found  from  Eq.   13: 


tan  (2(p  )= 


2   2 
1  ^  (1-Y2)  -  (X2-Y2)2 
Y^    V        NT 


M 


-[2(YZ-1-M2)+  4j-  ]-  2(Y-1) 


(65) 


To  find  the  equation  of  a  contour  curve  for  constant  argument  in  the   Y'~X 
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(oo-j8)  plane  we  use  the  notation    T  =  tan  (2cp)  and  expand  the  Eq.  65  in 

2  2 

powers  of   X      so  that  a  quadratic  equation  in   X      is  obtained. 


C4'X4  +  C2'X2  +  C^  =  0  •  (66) 


The  coefficients    C ',  C '/  C  '    depend  only  on  Y,  M,  and  T  and  are 
given  by  the  following  expressions: 


2 ',.2 


2    .„2,  2    ^.irj    ,    .,2,     M   I.^JUJl    ,x2.,,4 


C(J-tV[2(Y-1-MZ)+  ~]  -4T2(Y  -1)Y2[2(Y  -1-M2)+  ^-  ]+4T2(Y  -1)Z+M4      (67) 


C2'=  -2Y2T2[2(Y2-1-M2)+  ^.]2+4T2(Y2-l)[2(Y2-l-M2)+  ^]-  ^(1-Y2)     (68) 


C4'=T2[2(Y2-1-mV^    ]2  +  -Sl  (69) 


From  Eq.  66  we  obtain  the  explicit  form  of  the  equation  for  the  contour 
curve  for  constant  argument. 


2C4  V   V   2C4'    '  C4' 

Again,    X    must  be  real  and  positive.    Some  representative  contour  curves 
given  by  the  above  equation  are  shown  on  Figures    7  and  9. 
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Fig,   6    Contour  curves  for  constant  values  of  the  modulus    o    of  the  function 
u  =  o  ZE   in  the  complex  region  of  the    w  -  0    plane  at  the  value  of 
M  =  0.8. 
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Fig.  7    Contour  curves  for  constant  values  of  the  argument    <p     of  the 
function    U  =    o  ^   in  the  complex  region  of  the    w  -  £   plane 
at  the  value  of    M=0,8o 
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Contour  curves  for  constant  values  of  the  modulus  o  of  the 
function  U  =  £P  in  the  complex  region  of  the  b>  -  $  plane 
at  the  value  of    M  =  1  „  6 . 
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Fig.  9    Contour  curves  for  constant  values  of  the  argument    to   of  the 
function    U  =    o/!f_     in  the  complex  region  of  the    &j  -  # 
plane  at  the  value  of    M=l,6. 
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It  is  interesting  to  note  the  curves  delineating  various  zones  in 
the  complex  region.    On  tie  left  boundary  of  the  complex  region  between 
the  origin  and  the  upper  tangent  point  the  value  of  the  argument  of    U    is 
tt/2,  in  the  upper  tangent  point  the  value  of  the  argument  is  not  defined 
and  can  assume  any  value,  and  between  the  upper  tangent  point  and  the 
point  (X=l,  Y=l)  the  value  of  the  argument  is  zero.     In  this  point,  again, 
the  argument  is  not  defined. 

On  the  right  boundary  of  the  complex  region  between  the  origin 
and  the  lower  tangent  point  the  value  of  the  argument  is  zero,  in  the 
lower 'tangent  point  its  value  is  not  defined  and  between  the  lower  tangent 

point  and  the  point  (X=l,  Y=l)  the  value  of  the  argument  is    tt/2. 

2 

On  curves  where  Re(U  )  =  0  the  value  of  the  argument  is    7r/4. 

Contour  curves  converge  in  the  upper  and  lower  tangent  points,  in  the 
point  (X=l,  Y=l),  and  in  the  origin. 

It  is  now  possible  to  obtain  a  picture  of  the  behavior  of  the 

2  2 

functions    U,    and    U      in  various  regions  of  the    Y-X    (to  -  jS)  plane. 

The  curves  given  by  Equations  14,  15.  25,  26  and  the  lines  Y=l  and 


Y=  yl+M  divide  the    Y-X    (oj  -  j3)    plane  into  12  or  13  zones  depending 

whether  the  value  of    M    is  less  or  greater  than  unity.     It  has  already 

been  pointed  out  that  within  the  region  defined  by  Equations  14  and  15 

2  2 

the  functions    U,    and    U_    assume  complex  conjugate  values.     In  the 

2 

other  regions  the  following  has  been  found  concerning  the  sign  of    U, 

and    U  '. 
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Fig,   10    Zones  in  the    to  -  j8    plane  in  which  the  functions    U?    and    U? 


are  of  the  same  or  of  opposite  sign, 
values  of    M  <   1 „ 


General  behavior  for 
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Fig.   11    Zones  in  the    w-  S   plane  in  which  the  functions 

2  2 

Uj    and    U2    are  of  the  same  or  of  opposite  sign. 

General  behavior  for  values  of    M    >    1 . 
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At  Values  of    M    smaller  than  unity,  in  the  region  above  the  upper  branch 

2  2 

of  the  curve  given  by  Eq.  25.     Both    U1    and    LL    are  positive.     Between 

the  upper  branch  of  the  curve  given  by  Eq.   25,  and  the  curve  given  by 


/         2  2  2 

Eq.  26,  above  the  straight  line    Y=  y/l+M      ;    U      >  0,  U2    <  0.     On  the 

I         2-2  2 

line    Y=  \J  1+M         the  function    U?    has  a  pole  while  the  function    U. 

is  continuous.    To  the  right  of  the  curve  of  Eq.  26  and  above  the  line 

/         2  2  2 

Y=v/l+M         both    U,    and    U~    are  negative.     Between  the  lines 

I  9  2  2 

Y=l    and    Y=  J  1+M*     both    U,    and    U2    are  positive  to  the  left  of  the 

2 

curve  of  Eq.  26  while  to  the  right  of  the  curve  of  Eq.  26  the  value  of    U, 

'I 


2  2 

are  positive  to  tne  leit  ot  tne 

t  of  the  curve  of  Eq.  2  6  the  value  of    U 
is  negative  and  the  value  of    U9    positive.     Under  the  line    Y=l    and 


2 

above  the  curve  given  by  Eq.  26  the  value  of    U.     is  negative  while  that 

2 
of    XJ      is  positive.     Under  the  line    Y=l    and  under  the  curve  of  Eq.  26 

and  above  the  left  boundary  of  the  complex  region  given  by  Eq.   16  both 

2  2 

U,    and    U9    are  positive.     Under  the  curve  given  by  Eq.  26  and  to  the 

2 
left  of  the  boundary  of  this  complex  region  given  by  Eq.   16  both    U,    and 

2 
U_    are  negative.    Between  the  right  boundary  of  the  complex  region 

2 

given  by  Eq.   18  and  the  lower  branch  of  the  curve  of  Eq.  25  both    U,    and 

2 
U?    are  positive.     Under  the  lower  branch  of  the  curve  of  Eq.  25  the  function 

2  2 

U,     is  negative  while    U«    is  positive.      To  the  right  of  the  right  boundary 

of  the  complex  region  given  by  Eq.  18  under  the  line  Y=l    and  above  the 

2  2 

lower  branch  of  the  curve  of  Eq.  25  both    U,    and    U9    are  negative. 

2  2 

Figure  10  shows  the  zones  and  signs  of  the  functions    U,    and    U„    for 

the  case  when    M  <  1 . 
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For  values  of    M    greater  than  unity  the  lower  branch  of  the  curve  given 

by  Eq.  25  and  the  line    Y=l    intersect  and  modify  the  picutre  with  the 

2 

addition  of  another  zone.     In  suGh  cases  the  behavior  of  functions    U. 

2 

and    \J      is  modified  as  follows.    To  the  right  of  the  lower  branch  of  the 

2 
curve  of  Eq.  25  and  below  the  line    Y=l    the  function    U,     is  negative 

2 
while    U„    is  positive.     In  the  newly  formed  triangular  zone  bounded  by 

the  right  boundary  of  the  complex  region  given  by  Eq.   18,  the  line    Y=l , 

2  2 

and  the  lower  branch  of  the  curve  given  by  Eq.  25  both    U..    and    U^    are 

negative,  likewise  they  are  negative  in  the  region  above  the  line    Y=l 

and  below  the  lower  branch  of  the  curve  given  by  Eq.  25.     In  the  region 

to  the  right  of  the  curve  of  Eq.  26,  above  the  line    Y=l  and  above  the 


lower  branch  of  the  curve  of  Eq.  25  and  below  the  line    Y=  \f  M  +1 

2  2 

the  function    U,     is  negative  while  the  function    U~    is  positive.    The 

2  2 

behavior  of    U,    and    U_    in  the  other  zones  remains  the  same  as  for 

cases  when    M  <     1.    Figure  11  shows  the  zones  and  signs  of  functions 

2  2 

U      and    LL    for  this  case    (M>  1) . 

The  partition  of  the    03  -0     (Y-X)  plane  into  zones  as  described  above 
is  quite  general  and  is  valid  for  all  possible  cases.    These  cases  fall 
under  one  of  the  two  types,  viz.  ,  those  with    M<  1    and    those  with 
M>  1.    There  is,  however,  some  variation  in  the  sequence    of  special 
frequencies    (Y? ,  Y«  ,  — =-  +  \J  "~r~  +1'   vM  +1  '    etc«)  depending  in 
which  interval  lies  the  value  of    M.    The  sequential  order  of  these 
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special  frequencies  depend  only  on  the  values  of    M.    The  possible  cases 
are  listed  below. 

o  <  m  <ry-  (71) 


0  <  Y2  <  M    <  -  f  +  /^  +  1     <    Y3  <  1</m2+1    <f +/t +1  (72) 


=L  <   M   <  1  (73) 


/r 


°<\<-f+J4^  <y3<m<i  </m2+1  <¥+J4+i    ™ 


1  <  M    <  1.5  .  (75) 


0  <Y2<  -  f  +N/if  +1  <  Y3  <  1.  <  M  <v£T  <  f  +  7^  +1  (76) 


1.5  <  M  (77) 


o<-f +yMi+i<Y2<Y3<i<  M<yM2+i  <f  +/-T.*1     (78) 


2  2 

In  addition  to  the  general  behavior  of  the  functions    U.    and    UL    or, 

if  roots  are  taken,     U.    and    U_  ,  it  is  also  desirable  to  have  contour 
representations  of  these  functions  in  the    a)  -  j8    (Y-X)    plane  where    U. 
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and    TJL    have  either  real  or  imaginary  values,  i.e.  ,  outside  of  the  complex 
region. 

The  equation  of  contour  curves  in  the  to  -  0  (Y-X)  plane  for  which 
U  is  a  real  or  an  imaginary  constant  can  be  readily  obtained  by  solving 
Eq.  3  for  X. 


xW-u2*  m2u2+2y2(i-y2>  +    mV4y2(y2-i)(y2-i-u2)+m2u4 

2(Y2-1)  (Y2-M2)    ~  2(Y2-1)  (Y2-M2) 


This  is  a  simple  solution  of  a  biquadratic  equation    in   X.    For  contour 

curves  in  the  first  quadrant  of  the    w  -  fi    (Y-X)  plane  only  positive 

2 

values  of   X     are  acceptable.    There  appears  to  be  no  simple  rule  to 

determine  whether  the  curves  obtained  from  Eq.  79  are  the  contour 
curves  of    U,    or    U9 .    This  question  can  be  resolved  by  calculating  with 
the  obtained  values  of   Y   and    X    the  desired  function    U,    or    U?  ,  as  the 
case  may  be,  using  the  rule  concerning  the  selection  of  the  sign  as 
stated  above  and  then  determining  whether  the  calculated  value  is  equal 
to  the  value  of    U    originally  set  in  Eq.  79.     If  the  calculated  value  of 
U,     or    \J      is  equal  to  the  value  of    U    set  in  Eq.79  the  point  in  question 
lies  on  the  curve    U1  =  constant  or    U9  =  constant  respectively. 

Equation  3   can  also  be  solved  with  respect  to   Y   albeit  not  so  readily 
since  the  arrangement  in  powers  of   Y    leads  to  a  biquartic  equation. 
This  biquartic  equation  can  have  complex  roots  and  the  selection  of  proper 
curves  is  even  more  laborious  in  this  case.    There  appears  to  be  no  point 
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in  solving  the  Eq.  3  for  Y,  therefore,  since  the  solution  with  respect    to 
X    can  be  obtained  by  simple  algebra. 

Some  representative  contour  curves  for  various  real  and  imaginary 
values  of    U.    and    U-    are  shown  on  Figures  12,  13,  14,  and  15. 
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Fig,   12    Contour  curves  for  constant  real  and  imaginary  values  of    U. 
in  the    to  -  j3   plane  at  the  value  of    M  =  0.8. 
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Fig=   J 3    Contour  curves  for  constant  real  ana  imaginary  values  of    U 
in  the    ^  -  B   plane  at  the  value  of    M  =  0,8,  ' 
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cup 

Fig.   14      Contour  curves  for  constant  real  and  imaginary  values  of    TT 
in  the    us  -  B   plane  at  the  value  of    M  =  1„6„ 
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Fig.   15    Contour  curves  for  constant  real  and  imaginary  values  of    U2 
in  the    u>  -  B    plane  at  the  value  of    M  =  1.6. 
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